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Abstract 

We present a program that allows for the computation of tensor products of irre- 
ducible representations of Lie algebras A — G based on the explicit construction of 
weight states. This straightforward approach (which is slower and more memory- 
consumptive than the standard methods to just calculate dimensions of the tensor 
product decomposition) produces Clebsch-Gordan coefficients that are of interest 
for instance in discussing symmetry breaking in model building for grand unified 
theories. For that purpose, multiple tensor products have been implemented as well 
as means for analyzing the resulting effective operators in particle physics. 

Key words: tensor product; tensor product decomposition; Clebsch-Gordan 
coefficients; Lie algebra; multiple tensor product; model building; symmetry 
breaking; GUT. 
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Programming language: O'Caml 
Computer: i386-i686, x86_64 

Operating system: cross-platform, for definiteness though we assume some UNIX 
system. 

RAM: > 4 GB commendable, though in general memory requirements depend on 
the size of the Lie algebras and the representations involved. 

Keywords: tensor product; tensor product decomposition; Clebsch-Gordan coeffi- 
cients; Lie algebra; multiple tensor product; model building; symmetry breaking; 
GUT. 

Classification: 4.2, 11.1 

External routines/libraries: nums . cma (exact integer arithmetic, provided by the 
' Caml package) 
Subprograms used: none 
Nature of problem: 

Clebsch-Gordan coefficients are widely used in physics. This program has been writ- 
ten as a means to analyze symmetry breaking in the context of grand unified theories 
in particle physics. As an example, wc computed the singlets appearing in higher- 
dimensional operators 27(8)27(8)27(8)78 and 27(827(827(8650 in an E'e-symmetric 
GUT. 

Solution method: 

In contrast to very efficient algorithms that also produce tensor product decompo- 
sitions (as far as outer multiplicities / Littlewood-Richardson coefficients are con- 
cerned) we proceed straightforwardly by constructing all the weight states, i.e. the 
Clebsch-Gordan coefficients. This obviously comes at the expense of high memory 
and CPU-time demands. Applying Dynkin arithmetic in weight space, the algo- 
rithm is an extension of the one for the addition of angular momenta in su{2) ~ Ai, 
for reference see [1]. Note that, in general, Clebsch-Gordan coefficients are basis- 
dependent and therefore need to be understood with respect to the chosen basis. 
However, singlets appearing in (multiple) tensor products are less basis-dependent. 
Restrictions: 

Generically, only tensor products of non-degenerate or adjoint representations can 
be computed. However, the irreps appearing therein can subsequently be used as 
new input irreps for further tensor product decomposition so in principle there is 
no restriction on the irreps the tensor product is taken of. In practice, though, it is 
by the very nature of the explicit algorithm that input is restricted by memory and 
CPU runtime requirements. 
Unusual features: 

Analytic computation instead of float numerics. 
Additional comments: 

The program can be used in "notebook style" using a suitable O'Caml toplevel. 
Alternatively, an ' Caml input file can be compiled which effects in processing that 
is approximately a factor of five faster. The latter mode is commendable when large 
irreps need to be constructed. 
Running time: 

Varies depending on the input from parts of seconds to weeks for very large repre- 
sentations (because of memory exhaustion). 

References: 
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1 Introduction 



We aim at decomposing tensor products of irreducible representations (irreps) 
of simple Lie algebras A to G. By means of explicitly constructing weight 
states via successive lowering operations and complementing dominant weight 
spaces we compute Clebsch-Gordan coefficients. 

This work is organized as follows: in Section 2 we introduce our notations for 
Lie algebras and discuss their basic properties. Then, wc develop the machin- 
ery that is needed to explicitly decompose tensor products. In particular, we 
specify the basis vectors for which our algorithm yields Clebsch-Gordan co- 
efficients. Section 3 is an overview of the software structure whose individual 
components are described in Section 4. Installation instructions and test run 
examples are given in Sections 5 and 6, respectively. 



2 Theoretical background 



2.1 Lie algebras: basic facts and representation theory 



In this subsection we summarize the basic facts about the representation the- 
ory of Lie algebras to make the paper self-contained and to fix our notation. 
Most properties discussed here can be found in the hterature, e.g. [1,2,3,4]. 



A Lie algebra is usually defined via the fundamental Lie bracket. 



where the Lie algebra is completely determined by knowledge of its structure 
constants, Z"*"^. The Lie bracket fulfils the Jacobi identity. 



T", T^T^ + T\[T\T 



+ 



rj-}C rpa rpb 







(2) 



A Lie algebra is called simple if it is neither Abelian, i.e. if [g, g] ^ 0, nor has 
any non-trivial ideal, i.e. a subalgcbra t) C g with [i),g\ = f) other than the 
zero element and the algebra g itself. The simple Lie algebras can be classified 
by simple geometric methods [5,6], and can be grouped into the four infinite 
families of rank-n algebras, namely the special unitary algebras A„ = su{n + 
1), the special orthogonal algebras of odd dimension, 5„ = so{2n + 1), the 
symplectic algebras C„ = 5p(2n), as well as the special orthogonal algebras of 
even dimension, Dn = so{2n). Furthermore, there are five so-called exceptional 
Lie algebras, specified by their corresponding ranks, Eq, E^, E'g, F4, and 6*2. 



4 



A representation (rep) of a Lie algebra is a homomorphism from the set of 
generators to a corresponding set of linear mappings, Ry-, over a vector space, 
V, denoted in a compact way: {V^Rv)- The dimension of V is called the 
dimension of the representation. 

A maximal subset of mutually commuting Hermitian generators. Hi, i = 
1, . . . ,n is called a Cartan subalgebra. Its dimension is the rank of the Lie 
algebra. Its elements, the Cartan generators, fulfil 

Hi^Hj [Hi,Hj]^0. (3) 

Hence, they can be simultaneously diagonalized. As they form a linear space, 
we can choose a basis where tr [HiHj] = koSij, where the Dynkin index ko 
depends on the representation and the normalization of the generators. Any 
state in a representation D can be (up to possible degeneracies from other 
symmetries) uniquely labelled by the eigenvalues of the corresponding Cartan 
generators: 

Hi\fi,D) = fii\fi,D). (4) 

The eigenvalues //j are called weights, they can be grouped together in a vector 
of the dimension of the rank, called a weight vector. A scalar product on weight 
vectors a ■ ii = ai/ii is defined via index summation over i — 1, . . . , rank(0). 

A special representation is the adjoint representation, where the generators act 
by means of the Lie bracket on themselves, |T^) := | T"-,T^ ). Hence, the 
adjoint representation has the dimension of the rank of the algebra. Most of 
the facts about representations, weights etc. can be found in [2,3]. The weights 
of the adjoint representation are called roots. The non-Cartan generators, 
labelled by their root vectors a as i^^,, come in pairs of adjoints to each other 
with E_a — E'l- These form SU{2) subalgebras of the Lie algebra under 
consideration with Ea and E^a being the corresponding raising and lowering 
operators: 

[E^,E_^] = a-H . (5) 

This choice of operators is called the Cartan- Weyl basis. As the finite-dimen- 
sional unitary representations of SU{2) are well-known, all reps of all Lie 
algebras can be constructed. Let p be the maximum number a raising operator 
can be applied to an arbitrary state of a rep D, q the maximum number 

one can descend, then SU{2) algebra allows to relate the differences of the 
numbers {p — q) and {p' — q') for two different roots corresponding to two 
different raising operators to an angle between the roots, and it can be shown 
that this angle can only be 90, 120, 135, or 150 degrees. 

To specify a highest weight of a rep (like m = +j for spin j in quantum 
mechanics) one introduces an order on the space of weights by defining a 
weight to be positive if its first non-vanishing component is positive. With 
this help one can find a basis in the root space of a Lie algebra by taking the 
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so-called simple roots. These are those positive roots which cannot be written 
as sums of other positive roots. There are as many simple roots as the rank 
of the algebra. 



This positivity definition allows immediately for the ordering, namely two 
weights are ordered, /i > u ii /j, — u is positive. In the adjoint representation, 
the positive roots correspond to raising operators and the negative roots to 
lowering operators. The highest weight by definition cannot be raised by any 
raising operator. 

There is a diagrammatic notion of the system of simple roots, the Dynkin 
diagrams. Here, simple roots are denoted by circles, which are unconnected if 
the simple roots are orthogonal to each other, and connected simply, doubly 
or triply, if the enclosed angle between the two roots is 120, 135, or 150 de- 
grees, respectively Our conventions for the labelling of the roots are shown 
in Fig. A.l in Appendix A.l. Actually, there are extended Dynkin diagrams 
shown for all simple Lie algebras, augmented by the lowest root which is im- 
portant for the determination of the maximal subalgebras. Shorter roots (by 
a factor ^/2 for double hnks, a factor ^/3 for triple hues) are denoted by filled 
circles. 

The Cartan matrix allows to calculate the difference between lowering and 
raising possibilities for states corresponding to a positive root, (f) = J2j kja^: 



- = ^ kjAji with 
j 



Note that this is the definition in [2,3], but the transpose of [4]. The jth row 
of the Cartan matrix consists of the — values of the simple root . The 
Cartan matrices in our notation can be found in Appendix A. 2. 

By this means, it is straightforward to construct all possible reps of all simple 
Lie algebras, and there are many examples in the literature [2,3]. 

Any irrep is uniquely determined by its so-called Dynkin coefficients, £^ , which 
are defined via 

= . (7) 

Every set of & gives a /i which is the highest weight of one irrep. Hence, we 
choose as input a sequence of rank(0) non-negative integers to specify an irrep. 
Defining as fundamental weights those weight vectors for which 2a^ ■ fi'^/a^'^ = 
Sjk, every highest weight can then be uniquely written as = J2j^=i^^^ ^'^ 1^'' ■ 
The reps corresponding to those fundamental weights are the fundamental 
reps, again their number is rank(0). 




(6) 
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Fig. 1. Example for the construction of a rep of a Lie algebra, the fundamental 27 
of Eq, with Dynkin coefficients 100000. 

As an example we show the construction of all states of the fundamental irreps 
27 of Eq. We start from the fundamental weight 100000 as highest weight, and 
descend with the help of the Cartan matrix from it until we reach the state 
with the lowest weight at the bottom. What is shown in addition is the Dynkin 
coefficient of the lowest root. 

The program CleGo has a subprogram which is independently steerable by 
the user called lie . {opt I bin}. From the command line, the user specifies the 
Lie algebra under consideration as well as the irrep to be constructed. E.g. to 
produce the adjoint rep of SU{3) one uses ./lie. opt -su 3 -rep 11 which 
yields 



Lie algebra 




SU(3) 




Highest weight: 


(1,1) 




Dim. of irrep 




8 




1, L6v:0, Deg 


1 


(1,1), -2 


(0,0) 


2, Lev:l, Deg 


1 


(2,-l),-l 


(0,1) 


3, Lev:l, Deg 


1 


(-1,2),-1 


(1,0) 


4, Lev: 2, Deg 


2 


(0,0) ,0 


(1,1) 


6, Lev: 3, Deg 


1 


(1,-2),1 


(1,2) 


7, Lev: 3, Deg 


1 


(-2,1),1 


(2,1) 


8, Lev: 4, Deg 


1 


(-1,-1), 2 


(2,2) 



To reproduce the example from Fig. 1 of the irrep 27 of Eq one types in 
./lie. opt -e6 -rep 100000. This produces: 
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Lie algebra : E6 



Highest weight: (1,0,0,0,0,0) 
Dim. of irrep : 27 



1, Lev:0, Deg:l (1,0,0,0,0,0),-! (0,0,0,0,0,0) 

2, Lev:l, Deg:l (-1 , 1 ,0,0,0,0) , -1 (1,0,0,0,0,0) 

3, Lev:2, Deg:l (0,-1 , 1 ,0,0,0) ,-1 (1,1,0,0,0,0) 

4, Lev:3, Deg:l (0,0,-1 , 1 ,0, 1) ,-1 (1,1,1,0,0,0) 

5, Lev:4, Deg:l (0,0,0,1,0,-1) ,0 (1,1,1,0,0,1) 

6, Lev:4, Deg:l (0,0,0,-1 , 1 , 1) ,-1 (1,1,1,1,0,0) 

7, Lev:5, Deg:l (0,0, 1,-1 , 1 ,-1) ,0 (1,1,1,1,0,1) 

8, Lev:5, Deg:l (0,0,0,0,-1 , 1) ,-1 (1,1,1,1,1,0) 

9, Lev:6, Deg:l (0,0,1,0,-1,-1) ,0 (1,1,1,1,1,1) 

10, Lev:6, Deg:l (0, 1 ,-1 ,0, 1 ,0) ,0 (1,1,2,1,0,1) 

11, Lev:7, Deg:l (0, 1 , -1 , 1 ,-1 ,0) ,0 (1,1,2,1,1,1) 

12, Lev:7, Deg:l (1 , -1 , , , 1 , 0) ,0 (1,2,2,1,0,1) 

13, Lev:8, Deg:l (0, 1 ,0, -1 , , 0) ,0 (1,1,2,2,1,1) 

14, Lev:8, Deg:l (1 , -1 , , 1 , -1 ,0) , (1,2,2,1,1,1) 

15, Lev:8, Deg:l (-1 ,0 ,0 ,0 , 1 , 0) ,0 (2,2,2,1,0,1) 

16, Lev:9, Deg:l (1 , -1 , 1 ,-1 ,0,0) , (1,2,2,2,1,1) 

17, Lev:9, Deg:l (-1 , , , 1 , -1 ,0) ,0 (2,2,2,1,1,1) 

18, Lev:10, Deg:l (1 , ,-1 ,0,0, 1) ,0 (1,2,3,2,1,1) 

19, Lev:10, Deg:l (-1 ,0, 1 , -1 ,0 ,0) ,0 (2,2,2,2,1,1) 

20, Lev:ll, Deg:l (1 , ,0 ,0 ,0 ,-1) , 1 (1,2,3,2,1,2) 

21, Lev:ll, Deg:l (-1 , 1 , -1 , , , 1) ,0 (2,2,3,2,1,1) 

22, Lev:12, Deg:l (-1,1,0,0,0,-1),! (2,2,3,2,1,2) 

23, Lev:!2, Deg:! (0 ,-! ,0,0,0, 1) ,0 (2,3,3,2,1,1) 

24, Lev:13, Deg:l (0,-1 , 1 ,0,0,-1) , 1 (2,3,3,2,1,2) 

25, Lev:14, Deg:l (0,0,-1 , 1 ,0,0) , 1 (2,3,4,2,1,2) 

26, Lev:15, Deg:l (0,0,0,-1 , 1 ,0) , 1 (2,3,4,3,1,2) 

27, Lev:16, Deg:l (0,0,0,0,-1,0) ,1 (2,3,4,3,2,2) 



The program repeats the Lie algebra, the highest weight chosen, calculates 
the dimension of the irrep, and then constructs it. In each line, the entries are 
the number of the state, the level below the highest weight (with the highest 
weight having level zero), the level of degeneracy (here there is only the trivial 
degeneracy of the rank(0) zero weights of the adjoint). The two rightmost 
columns show the Dynkin coefficients of the corresponding weight with the 
Dynkin coefficient of the lowest root split by a comma, and the number of 
applications of the corresponding lowering operators from the highest weight. 
E.g. state number 4 of the 27 of Eq is generated by applying the raising 
operators corresponding to the simple roots a^'^'^ exactly once. 

While the construction of the weights of an irrep is just simple addition of rows 
of the corresponding Cartan matrix to the highest weight, the determination of 
the dimension of an irrep can be done by means of the system of the positive 
roots of the Lie algebra q, A_|_. The roots of a Lie algebra can be uniquely 
decomposed (up to Weyl symmetries of the root lattice) into the zero roots, 
the positive and the negative roots. Our conventions for the positive roots 
can be found in Appendix A. 3, again in the Dynkin basis. In the following 
discussion we denote the highest weight of an irrep in the Dynkin basis by 



The dimension of an irrep R can be derived from characters and is given by 
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Weyl's formula [7]: 



dim(i?) 



n 

aeA+ 



(A + S) - a 
6 ■ a 



n 



raiik(g) 



= 1 



^rank(g) , 



+ 1 



(8) 



where S — ^ Z^aeA+ ck- ^? ^-re the coefficients of the root a with respect to 



the simple roots, a = Y.^^^=i^^^ kfa^, rij are the Dynkin labels of the highest 

weight, A = "j/^^- and coj is the weight of the simple root in the 

Dynkin diagram. The weights of the roots are their squares, they are one for 
all roots in Lie algebras with only simple links, one for the shorter and two 
for the longer roots in S0{2n + 1), Sp{2n), and F4, and one for the shorter, 
three for the longer root in 6*2. 

The degree of degeneracy for a specific state in an irrep can be calculated 
by the Freudenthal recursion formula [8] recursively from the degeneracies of 
all states of higher levels than the corresponding state. The formula goes as 
follows: for an irrep R with highest weight A the multiplicity m(A) of a possible 
weight A = A — QjOi-' with qi,q2, . . . , ?rank(0) all non- negative integers is 

given by: 



[(A + 5) ■ (A + 5) - (A + 5) ■ (A + 5)]m(A) = 2 ^ ^ m(A+fca)((A+A;a) -a) 

aeA+ k 

(9) 

The second sum on the right hand side extends over all those values of k for 
which X + ka whose level is less than that of A. 5 = | J2aeA+ a as in the Weyl 
formula. The Freudenthal formula can be recast in a more usable form which 
has been used in the implementation in our program: 



' rank(g) 

i=i 



Wo 



1 



rank(g) 



no 



i=l 



m(A) 



{rank(g) 



rank(g) 

nj + E (kkt - q,) 
1=1 



(10) 



Here, A 



E;=f^^n,/.^ 



A 



^_^rank(g) ^ 



a 



^rank(g) j^^^j ^ 



■jj=L '"jr' 1 " ^j=l 

weight of the simple root a' in the Dynkin diagram. The second sum is defined 
in the same way as the one in (9). The highest weight has m(A) = 1 which 
allows to obtain the degeneracies or multiplicities first for level 1, then for 
level 2, and so on. For example, every level-1 weight has the form A = A — a-', 
where is a simple root. The only non- vanishing term on the right-hand side 
of (10) is for a — and k — 1, which is 2m(A)A ■ — A • aK In the next 
step, the degeneracies of the level- 2 weights are determined from those of level 
and 1, and so on. 



Regarding the determination of the Dynkin coefficient of the lowest root, 
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we have to determine first the Uncar combination with respect to the 
roots. These are (in our conventions): 



simple 



Bn 
Cn 
Dn 

Eg 
G2 



a — 



a 



a 



a 



a 



a 



a 



a 



a 



i=l 

n 
1=2 

n-1 

-2 ^ a* - a" 



-a 



1=1 
1 



n—3 

2E 



a 



a 



n-l 



a 



-a^ - 2a^ - - 2a^ - - 2a^ 
-2a^ - 3a^ - 4a^ 
-2a^ - Aa^ 
-2a^ - 
-3a^ - 2a^ 



- - Aa^ 
2,a^ - 2a^ 



- 2a' - 3a^ 



(11a) 
(lib) 
(11c) 

(lid) 

(lie) 
(llf) 

(llg) 
(llh) 

flli) 



To calculate the Dynkin coefficients according to the formula 

■ n 



a 



o2 



(12) 



the formulae (11) are the same for the Dynkin coefficients for the series An 
and Dn as well as for -^6,7,8- For F4, and are shorter by a factor \/2 
as the other three, leaving us with the relation = — £^ — 2^'^ — 3£'^ — 2£^. 
In the case of G2, oP and are longer by a factor -\/3 as a^, which yields 
= — 2£^. For the two infinite series the last coefficient has to be 
rescaled by two, while for C„ all coefficients except for the last one have to be 
rescaled by two. In general, is always among the longer roots. 



2.2 Tensor product representation 



The definition of the tensor product of representations of a Lie algebra is 
as follows: Let be a Lie algebra and (V,i?i/); iW-iRw) representations of 
Q. We then construct another representation, the so-called tensor product 
representation, (V ® W, Ry <E) Rw) via 



{Rv ® Rw)ix) -.V 0W 
{{Rv <S> Rw){x)) {v (g) w) := {Rv{x)v) <S>w + v<^ {Rw{x)w) 



(13) 
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for all X G g. li the original representations are finite- dimensional and irre- 
ducible, the tensor product representation is also finite-dimensional which in 
the case of (semi)simple algebras implies that it is completely reducible 

n n 

i=i 1=1 
with irreducible representations {Ui,Rui)- 

The tensor product is commutative up to isomorphism, i.e. there is an isomor- 
phism f: V'^W-^W'^V such that 

Rw^v{x) o f ^ f o Rv^w{x), Vx e 0. 

Also, multiple tensor products of representations are associative up to isomor- 
phism which again means that there is an isomorphism 

f :U ®{V ®W) ^ {U ®V)®W 

that intertwines representations {U ® {V ® W), Ru <8) {Ry ® Rw)) and {{U (g) 
V) ® W, {Ru ® Rv) ® Rw) of 0. 

In particular, for a simple Lie algebra and finite-dimensional representations, 
we have the decomposition of the tensor product module 

V^a®^a' = 04a'^a. (14) 

i 

where modules V\ are uniquely specified by their highest weight A, etc.. The 
number of times an irrep appears on the right-hand side of (14) is given by 
which are non-negative integers to be denoted tensor product multiplicities or 
Littlewood- Richardson coefficients of g. 

The Cartan-Weyl basis of the Lie algebra is translated to operators in the 
respective representations by means of the algebra homomorphism of the rep- 
resentations. Let its Cartan generators be represented as {Hvjji and {Hv^,)i 
on Va and Va' respectively, the definition of the tensor product representation 
in (13) shows that the Cartan generators acting on the product space read 

H,= {Hy^)i®l + l®{Hv^,)i. (15) 

Likewise, raising and lowering operators {Evj^)±a and {Ey^,)±a are given by 

E±a = (£;y)±a ® 1 + 1 ® {Ew)±a (16) 

where a denote the roots of the Lie algebra g. One can easily check that the 
generators of the product representation again form a Cartan-Weyl basis. In 
particular, we have 
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[Hi, E±a\ — icKjE'ia. 



(17) 



It is therefore again appropriate to discuss the tensor product representation 
in terms of weights and roots, that is, in Dynkin weight space. Prom (15), it 
is clear that all weights of the tensor product representation are obtained by 
adding weights of the respective irreps the tensor product is taken of. The 
multiplicities of these weights A is 

multy^^v-A/l-^) = J2 multvx(/x) • multv^,(/i')- 

with ii+ij,'=\ 

It then follows that the highest weight in V\ ® V\i is given by the sum of the 
highest weights of Va. and Va.', i.e. A + A', and appears with multiplicity one. 



2.3 Explicit lowering 

Before describing the decomposition algorithm we need to discuss in some 
detail how the Cartan-Weyl generators act on representation space. This dis- 
cussion is close to the one in [9]. Given an irrep of a simple Lie algebra we 
choose basis states in the weight subspaces that, in general, are degenerate. 
In what follows we use the following basis states 

\{w + a)v), where T = 1, . . . , D^+a, 

l^c), where c = 1, . . . , D^j, 
\{w - a)c), where C = 1, . . . , D^^^, 

where denotes the degeneracy of the subspace of weight w. Those basis 
states are in general non-orthogonal. Without loss of generality, we want them 
to be always normalized to unity, i.e. 

{Wc\Wc) = 1, Vc. 

For every weight w in the representation, let M'^ be a matrix with Af^ = 
{wa\wb) and the matrix with G^f, = ((M"')~^)ab, then the identity operator 
in the degenerate subspace is given by 

li^) ^Y.Gl\wa){w,\. (18) 

a,b 

It helps to consider states 



12 



a 
b 

which satisfy 

{Wc\Wa) = (WalWc) = Sac- (19) 

Generically, raising or lowering of a degenerate weight state is a hnear combi- 
nation of degenerate states of the adjacent weights: 

r 

E-a^Wc) = J2 N-ai,Wc^{w-ai)c I - ai)c) (20) 

c 

Using (19) we can single out the lowering coefficient 

N_oci,Wa^{w-ai)A = IZG'ab"M(^^ - ai)B\E_ai\Wa)- (21) 
B 

From Ea = E'L^ and the defining relation in (20) we obtain 

A,B 

= {Wb\EaE_a\Wa) = (Wf,| + E_aEa\Wa) 

= {wb\wa){a,w) + J2GrA" {^b\E-a\{w + a)r){iw + a)A|^aka) 
r,A 

r,A,c,d 

where we made use of [E^, E_a] = (ct, H) and the identity operator from (18) 
in the next to last step, as well as (21) in the last step. Note that in the Dynkin 
basis, {ai,w) denotes the i-th Dynkin coordinate of w, i.e. Wi. 

Hence, setting a = 6 we obtain 

A,B 

= {a,w) + 

(22) 



r,A,c,d 



which is close to a recursion relation for the lowering normalizations: Given 
that one knows how to lower states of weight w + a and one also has the scalar 
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products of the basis states of both weight w and w + a, the right-hand side of 
(22) is completely determined. We see that normalizations for lowering states 
of weight w and scalar products of states of weight w — a need to be chosen 
consistently, that is, in such a way that (22) holds. 

There are two special cases in which the consistency equation in (22) simplifies 
considerably. First, in the case where the irrep is non-degenerate, we have by 
construction 

rw /„,, L,, \ 1 /^w I ( l\iw\ — l 



MZ = {wa\w,) = 1, G2 = ((^")"'J„, = 1, = 
which simplifies (22) to 

I = {a,w) + \N_a,{w+a)^w\ ■ (23) 

Thus, for non-degenerate irreps lowering normalizations are recursively de- 
termined up to the choice of phase factors. It is convenient to always take 
lowering phase factors to be unity. If the right-hand side turns out to be zero 
or negative, the state \w — a) does not exist. 

Secondly, we consider the case where the irrep is the adjoint representation 
which has a degenerate zero weight space but is otherwise non-degenerate. 
Starting with lowering the highest weight, normalization constants are the 
same as in the non-degenerate case until we are to lower the simple roots. For 
a simple rootQ w = lowering can only occur via a = a* and the consistency 
equation becomes 

rank 
A,B=1 

Now the basis in the subspace of weight zero is chosen in such a way that basis 
states therein result from lowering simple roots and are normalized to unity, 
i.e. 

. rv* rv*— ^0 1 if ^ = Z 

iV_,.,,.^o, = ' " (24) 



0, else 



and obtain 



^-a\a»-i-oJ — 2. (25) 



^ Upper indices label distinct simple roots in order to distinguish them from lower 
indices that label degeneracies. 
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Note that in (25) \N_ai^2a^^aA^ vanishes because 2a^ is never a root. Now, if 
w — and a — a^, we have 



= 2(0„|0,)2 = ^i^ (26) 
where (24) and (25) have been used and the last step is due to 

(OJO,) = ^^^^ (27) 

with the Cartan matrix A. Relation (27) will be shown in 2.3.1. Now let 
w — —a" and a — a^, then the condition in (22) reads 



|-^— a**,— a"— 



2 _ 



{a^, -a") + " )^A-^-a^(a''-a»)^->-Q;°-^-a^(Q;''-a°)A->-a" 

r,A 

which for a = 6 becomes (using (26) and unity phase factors) 

|-^-a",-Q;'»-^-2a<' 1^ = ~2 + E(G'°)rA-^-a",Or^-a"-^-a°,OA-^-a'' 

r,A 

= -2 + E(G^°)rA(MO)r„(MO)Aa|A^-a»,a<'^o„r 
r,A 

= -2 + (M°)aa|iV_a«,a.^0.r = -2 + 2 = 

while for a 7^ 6 it is simply 

I A^_Q,6^_Qa_^_^o_„6 1^ = (a'', — Q;") + 

since the negative simple root, —a", can only be obtained from lowering zero 
weights by the a-th root. Further lowerings are analogous to the ones in non- 
degenerate irreps. Combining all those cases one finds a nice recursion relation 
for lowering normalizations of the adjoint representation: 

Again, for simplicity, phase factors are all taken to be unity. 

Until now, we only specified lowering normalizations in non-degenerate or ad- 
joint representations. It is hard to find a general lowering scheme but, fortu- 
nately, this is not needed as one can built up arbitrary irreps from (multiple) 
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tensor products of non-degcnciate and adjoint representations and transfer 
bases and normalization constants. 



2.3.1 Scalar products of basis states of adjoint zero weights 

We now return to the scalar product (27) of the basis states of weight zero in 
the adjoint representation. Prom (20) and its Hermitian conjugate, we find 



\a 



= ^ {a^\ {[E^i,E_^i\ + E_^iE^i) \a') (29) 

where also (25) has been used. As — a* is an adjoint weight only if i = j, 
the commutator [E^j, E_a,i\ is non- vanishing only in the case i — j and is then 
given by 

[E^.,E_^^ = {a\H). 

For i = J, Eai\a'^) is zero because 2q;' is never a root. Hence, in this case we 
find 

(0,|0,) = ^ {a'\{a\H)\a') = \{a'\a') = 1 (30) 
On the other hand, for i j, we are left with 

{Oj\Oi) = ^{a^E_^iE^,\a'). (31) 

We already know that in this case we have 



E_ai \-a^) = ^ {a\ -a^) \- -a^) = ^J\Aji\ | - a* - a^) 
which can be used to fix the following raising normalizations: 



E^.\a^) = ^\A,,\\a' + a^) (32) 
It is now by means of (32) and its Hermitian conjugate that we end up with 



(0,|0,) = i(a^|^_,.^,.|a^) 



= + ^.|^. + ^ (33) 

Eventually, the two cases in (30) and (33) can be put together to what we 
wanted to prove, (27): 

/n In \ V A-abA-ba 
{(Ja\(Jb) - 7, 
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Note that this relation generahzes (Oa|0;,) = |Aa5|/2 which is the one from 
Anderson & Blazek [9]. It is apparent that their relation is restricted to the 
case where the Cartan matrix is symmetric and, thus, applies only to Lie 
algebra classes A, D and E. In contrast, our relation applies to all Lie algebras 
A to G. 



2.4 Decomposition algorithm 

There are various algorithms for the computation of tensor product decompo- 
sitions of simple Lie algebra representations. If it were only for the Littlewood- 
Richardson coefficients there exists very efficient means. For instance, there 
is a fast algorithm based on characters and Klimyk's formula [13] H] However, 
not only are we interested in multiplicities but also in the Clebsch-Gordan 
coefficients, that is, the coefficients of tensor product states in terms of the 
states of the two irreps the tensor product is taken of. We therefore proceed 
straightforwardly by explicitly constructing states in the modules of the tensor 
product representation. This, basically, generalizes the Clebsch-Gordan algo- 
rithm for SU{2) which is well-known from the addition of angular momenta. 
Such an explicit algorithm is also found in the works by Koh et al. [10] and 
Anderson & Blazek ([11], [12], and [9]) in which tensor products of some lower- 
dimensional irreps of Eq were computed. It will become apparent that in this 
approach we will need knowledge of consistent lowering normalizations and 
bases in the irreps the tensor product is taken of. 

Before presenting the algorithm we still need another definition: Dominant 
weights are weights in the product representation with non-negative Dynkin 
coordinates only. They are of interest as they can serve as highest weights of 
irreps in the decomposition. (Apparently, highest weights of irreps are always 
dominant weights.) Now, our Clebsch-Gordan algorithm is as follows: 

(1) Compute the highest weight as the sum of the highest weights of the two 
irreps. 

(2) Descend one level, order by weights and drop states until the remaining 
ones become linearly independent, i.e. they now form a basis. This is the 
CPU-intensive part. 

(3) Descend next level . . . until the full irrep is constructed. 

(4) Find dominant weights and list them. 

(5) Compute basis of dominant weight states in the tensor product repre- 
sentation from the weights and their degeneracies in the two irreps the 
tensor product is taken of. 



^ An implementation of this algorithm is found in a software called LiE which 
decomposes nearly any tensor product in less than a second [14] . 
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(6) Choose dominant weight with the most level^ and check if the comple- 
ment (of dominant weights from step 4 with respect to the ones from 
step 5) is non-empty. If so, this is the highest weight of the next irrep in 
the decomposition. Otherwise, remove this dominant weight from the hst 
and repeat step 6 as long as the list is non-empty. 

(7) Descend irrep whose highest weight is given by the dominant weight from 
step 6. Go to step 2. . . 

(8) Algorithm ends automatically, check dimensions. 

As an illustration of the algorithm, let us consider the following neat example 
for SU{3): 

3® 3 =? 

Both the 3 and 3 are non-degenerate and, thus, lowering normalizations are 
taken from (23) 

3: |10) ^ 1 ■ 111) 4 1 ■ |01) 
3: |01) 4 1 ■ 111) 4 1 ■ |10). 

Here, numbers over arrows denote the lowering operators corresponding to the 
respective simple roots. The highest weight in the tensor product representa- 
tion is given by 

|11) = |10)|01) 

from which we descend, using (16), and build up the essential part of what 
turns out to be an octet: 

|11)4|12) = |11)|01) A V2- IOO2) = |01)|01) + |ii)|ii) ^4 ... 

Ill) A |2i) = |io)|ii) 4 72 • |00i) = |io)|io) + |ii)|ii) '4... 

We need not descend further as dominant weights cannot occur below the level 
with the zero weight. Dominant weight states other than |11) are 

V2- |00i) = |10)|10) + |11)|11) 
72-1002) = |01)|01) + |11)|11). 

Now we list basis vectors for the |00)-space — |10)|10), |11)|11), |01)|01) — 
which turns out to be of dimension 3. Having already two states of this weight 
there can only be one additional state which is orthogonal to the first two. We 
find 

73 ■ |00) = |10)|10) - |11)|11) + |01)|01) 

^ For Lie algebras A to G there is an r-tuple, the so-called level vector, whose 
standard scalar product with the highest weight gives the level of the lowest weight, 
i.e. the maximal level [3]. 
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which, in fact, turns out to be a singlet. Note that all kets are normalized to 
unity. The result is well known, 

3(8)3 = 8 + 1. 

It is by the very nature of this algorithm that for larger Lie algebras (rank 
> 2) or representations with higher dimensions (> 10) computations grow 
rapidly in length and therefore become very tedious. The implementation of 
this algorithm is described in the following sections. In principle — that is, up 
to restrictions coming from memory and CPU-time demands — it enables us 
to compute any tensor product decomposition of irreps of Lie algebras A to 
G, as well as the corresponding Clebsch-Gordan coefficients. Note, however, 
that in general Clebsch-Gordan coefficients are not universal when arbitrary 
irreps with degenerate weight spaces are involved as they only make sense if, 
at the same time, bases in the degenerate weight spaces are specified. 
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3 Overview of the softwEire structure 



Our code is divided into the following parts: 



File 


Siumiiary 


Aux . ml 


Auxiliary functions such as list manipulation rou- 
tines (beyond the ones from O'Caml standard 
module List), combinatorical functions needed in 
Liealg.ml or simplification methods for big_int 
numbers. . . 


Liealg.ml 


Subprogram that is capable of computing weight 
lists and weight degeneracies 


Alg . ml 


Subprogram to model an algebra 


CGC.ml 


Main program for decomposing tensor product 
representations and computing Clebsch-Gordan 

coefficients 


CGC.nb.ml 


User notebook specifying input and output of the 
computation. 



Our software makes use of the modules List and Array from the O'Caml stan- 
dard library. Moreover, module Num from the Num library is used for arbitrary- 
precision integer / rational arithmetic. 



It is by the very nature of representation theory that the complexity of weight 
spaces in tensor product representations grow rapidly when larger irreps or 
irreps of larger algebras are involved in the tensor product representation. 
In such cases, it often does not make sense to output the resulting Clebsch- 
Gordan coefficients. Instead, one better proceeds until computations reach a 
much simpler result which may then be a suitable output. This is the reason 
why our program cannot be distributed as a single compiled user-friendly bi- 
nary executable. It should rather be understood as a package that provides 
the tools to do computations that require knowledge of Clebsch-Gordan co- 
efficients. We therefore included a notebook file CGC_nb . ml where such input 
and output can be specified. 
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4 Description of the individual softweire components 

4-1 Aux 

The module Aux. ml contains auxiliary functions and operations needed for 
the explicit construction of irreps and the Clebsch-Gordon decomposition of 
multiple tensor products. Parts of it bear similarities to the corresponding 
modules on lists in the O'Mega matrix element generator for quantum- field 
theoretical amplitudes [15]. This module contains e.g. functions to generate a 
list of all inequivalent pairs of numbers smaller than n. Most functions in Aux 
are rather primitive and may go without further explanation. 

We should, however, mention that here wc define the following tree structure, 
type 'a tree = Empty I Node of int * 'a * 'a tree * 'a tree, 
which will be used in the construction of multiple tensor products (see 
class irrep_in_tp in CGC.ml). Some tree manipulations are also contained 
in this module. 

4-2 Liealg 

This subprogram generates by explicit construction the list of all possible 
weights of any irrep in any simple Lie algebra. It determines all the states of an 
irrep either in the Dynkin basis or as a tuple of numbers how many times one 
has to descend from the (unique) highest weight to get to the corresponding 
weight of the irrep. A trivial piece of information is the level, which is just the 
number of total descends from the highest weight. The construction proceeds 
simply by consecutive subtraction of the rows of the Cartan matrix from the 
highest weight. In that construction, some states are constructed more than 
once (if they can be reached on different paths of descent from the highest 
weight); multiple states are eliminated and the states at each level will be 
sorted. 

The recursive Preudenthal formula (cf. (9), (10)) is used to determine the de- 
gree of degeneracy of each level, an information of importance for the explicit 
generation of the space of states of each irrep. 

Here, we briefiy describe the setup of the elementary objects and most im- 
portant functions of our software implementation in the module Liealg. A 

weight is denoted by a pair, where the first entry itself is a pair, giving the 
level in an irrep starting with zero for the highest weight, the second entry is 
a list of integers which has as nth entry the number of how many times the 
simple root a" has been applied to the highest weight. The second entry of 
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the weight pair is again an integer list, namely the Dynkin coefficients. The 
function complete_descent : liealgebra -> 

int list -> ((int * int list) * int list) list 
has as arguments a Lie algebra (being one of the elements A of int, B of 
int, C of int, D of int, E6, E7, E8, F4, or G2) and an integer hst giving the 
Dynkin coefficients of the highest weight of the representation. The result is 
the list of all weights of the irrep. 

This set of all weights is used by the function 

freudenthal: liealgebra -> int list -> (int * (int * int list) 
* (((int list * int) * int) * (int * int list)) list) list 

which has as input again the Lie algebra and a highest weight of an irrep, and 
gives back a list of all the weights (degeneracy, level and Dynkin coefficient), 
and a list of all states relevant for the recursive calculation according to the 
sum on the right-hand side of (10). It is a constructive recursive function 
implementing the Freudenthal algorithm. 

These two functionalities arc combined in the function: 
f ull_descent : liealgebra -> int list -> 

((int * int * int list) * int list * int) list 
which gives back just the list of weights for the irrep defined by the highest 
weight of type int list including the degeneracies. Each weight is augmented 
by the Dynkin number that corresponds to the zeroth root. 

When using the ' Caml toplevel a more readable output of the weight system 
is given by means of the function 

weights: Liealg. liealgebra -> int list -> unit 

where again the two input parameters correspond to the Lie algebra and the 
heighest weight. 

We just mention the function 

cartan: liealgebra -> int list list 

which gives the Cartan matrix for the corresponding Lie algebra. 



4.3 Alg 

The construction of explicit states in the tensor product representation re- 
quires some linear algebra techniques most of which are based on a suitable 
implementation of the Gaussian algorithm. However, as we aim at producing 
exact Clebsch-Gordan coefficients — as opposed to float numbers — we flrst 
have to model objects that represent numbers that are linear combinations of 
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square roots of non-negative integers with rational coefficients, i.e. 

oo 

E?"V^ (34) 

n=0 

where aU but finitely many Qn ^ Q are zero. We then extend the numbers in 
(34) to a field and also model the notion of a vector space over this field. Our 
linear algebra routines are then based on the latter two. 

We have the following classes: 

• class rational: 

Rational numbers Q based on the arbitrary-precision rational numbers from 
the Num library. Its constructor takes two arguments referring to the two sig- 
natures (numerator : int) and (denominator :int). Commonly used meth- 
ods are: 

add: rational -> rational 
multiply: rational -> rational 
divide_by: rational -> rational 
is_zero: -> bool 
is_bigger_than_zero : -> bool 
print : -> string 

• class linear_comb_rational_sqrts: 

Models algebra over the rational numbers generated by the square roots of 
non-negative integers, see (34). Its constructor takes an argument of type 
(rational * int) list where every list entry corresponds to a summand 
in (34). Important methods are: 
get_list: -> (rational * int) list 

add: linear_comb_rational_sqrts -> linear_comb_rational_sqrts 
multiply: linear_comb_rational_sqrts -> 

linear_comb_rational_sqrts 

remove_zeros : -> linear_comb_rational_sqrts: To remove list entries 
that correspond to zeros in (34) which can be remnants of methods add and 
multiply. 

symplify: -> linear_comb_rational_sqrts: Simplifying by extracting 
roots as much as possible and simplifying terms that include identical square 

roots. 

is_zero: -> bool 

print: -> string: String output to be used in the O'Caml toplevel. 

tex: -> string: Outputs number in T^format. 

export _math : -> string: Outputs number in Mathematica format. 

• class lcrs_field: 

Models the field extension of linear_comb_rational_sqrts numbers via 
the inclusion of fractions of linear_comb_rational_sqrts numbers. Its 



23 



constructor takes two linear_comb_rational_sqrts objects referring to 

the numerator and denominator. Methods of importance are: 

get_num: -> linear_comb_rational_sqrts: Yields the numerator of the 

fraction. 

get_denom: -> linear_comb_rational_sqrts: Yields the denominator of 
the fraction. 

remove_zeros : -> lcrs_f ield: Removes zeros in both the numerator and 
the denominator. 

symplify: -> lcrs_field: Symplifies the numerator and the denomina- 
tor. In the case where the denominator consists only of one term further 
simplifies the fraction, 
add: lcrs_field -> lcrs_field 
multiply: lcrs_field -> lcrs_field 
invert: -> lcrs_field: If non-zero, gives the inverse. 
is_zero: -> bool 
is_bigger_than_zero : -> bool 
print : -> string 
tex: -> string 
export _math: -> string 

• class ['a] lcrs_vector: 

Models the vector space over lcrs_f ield with basis vectors labeled by val- 
ues of type 'a. The constructor takes an argument of the following type: 
(lcrs_field * 'a) list where the list represents a linear combination of 
basis vectors of type 'a with lcrs_field coefficients. Useful methods of 
this class are: 

add: lcrs_field -> 'a -> 'a lcrs_vector: Adds another term to the 
list. 

add2: 'a lcrs_vector -> 'a lcrs_vector: Adds another vector of the 
same type ' a and returns the sum. 
get_list: -> (lcrs_field * 'a) list 
print: -> (string * 'a) list 

scale: -> lcrs_field -> 'a lcrs_vector: Scales the vector by a num- 
ber of our field. 

remove_zeros : -> 'a lcrs_vector 
symplify: -> 'a lcrs_vector 
is_zero: -> bool 



Based on lcrs_f ield, we now have the following functions: 

• number: int -> int -> int -> lcrs_field: 

Instantiates (via number a b n) objects of lcrs_f ield that represent num- 
bers of the simple form 

• gcd_of _lcrs_f ields : lcrs_field list -> lcrs_field: 
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Computes the greatest common divisor of a set of lcrs_f ield numbers. 

array_product : .. array -> .. array -> lcrs_field: 

This function takes two arguments that e.g. are of typ^ lcrs_f ield array 

and returns the scalar product provided that the two arrays have equal 

length. 

matrix_multiply : . . array array -> . . array array -> 
. . array array: 

Multiplies e.g. two matrices of type lcrs_f ield array array and yields a 
matrix of the same type, provided that dimensions of matrices fit. 



• gauss: lcrs_field array array -> lcrs_field array array 
-> lcrs_field array array: 

Gaussian algorithm (adapted from [16]) transforming a lcrs_f ield matrix 
to Gaussian shape. The second argument is meant to be a column vector 
representing the right-hand side of a system of linear equations. 

• solve: ... array array -> lcrs_field array: 

Gives an explicit arbitrary solution (provided there is one) to a system of 
linear equations which has already been brought into Gaussian shape. Here, 
the argument . . . (which represents the system of linear equations) is again 
to some extent polymorphic. It is sufficient to state that e.g. the output of 
gauss is of the right type. 

• invert: lcrs_field array array -> lcrs_field array array: 
Inverts (if possible) a square matrix over lcrs_f ield. Note that although 
invert is not needed in computing Clebsch-Gordan coefficients, it is useful 
when it comes to basis changes in CGC_nb. 

• linearly_dependent : .... list -> bool: 

Tests whether or not a set / list of vectors are linearly dependent. The 
argument .... is again polymorphic. As an example, an argument could 
be of type ['a] lcrs_vector list. Note that although this function is 
not needed in computing Clebsch-Gordan coefficients it proves useful in 
constructing bases etc. 



4.4 CGC 



This part contains the algorithm to explicitly decompose tensor products of 
irreducible representations of Lie algebras A - G. 

We have the following functions / variables: 

• level_vector : Liealg. liealgebra -> int list: 

^ This function is polymorphic in that the arguments are only restricted to be an 
array of an object that comes with suitable methods add and multiply. This is 
what we mean by . . above. 
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Returns the level vector of the Lie algebra. They are taken from [3]. As 
already stated in Section 2.4 dominant weights in the tensor product repre- 
sentation are sorted according to their maximal levels, the latter of which 
can be computed by means of the level vector. 

• adjoint: Liealg. liealgebra -> int list: 

Returns the highest weight of the adjoint representation for a given Lie 
algebra. Generically, we only have consistent lowering normalizations for 
non-degenerate or adjoint irreps, see Section 2.3. It is by means of this 
function, adjoint, that we can recognize adjoint representations!^ 

• empty.input = ( [] , (([], [] ) , [])): 

When instantiating an object of the class irrep (see below) the variable 
empty_input as an input parameter indicates that the standard lowering 
normalizations are used. This, of course, works in the case of non-degenerate 
or adjoint representations only. 

• scp_zero_weights : Liealg. liealgebra -> int -> int 
-> Alg. lcrs_f ield: 

Gives the scalar products of the basis states of the zero weight subspace in 
the adjoint representation from the Cartan matrices according to (27). Our 
basis states are defined in terms of descending the simple root states. The 
second and third argument of type int are restricted from 1 to the rank 
of the Lie algebra. The knowledge of those scalar products are needed for 
consistent lowering normalizations for adjoint representations, see (26). 

We have the following classes: 

• class ket: 

Models a generic state. Its constructor takes two arguments of which the 
first is an int list to specify the weight of the state and the second is of 
type int to label degeneracies. One of its methods are: 
print : -> string: Outputs a string containing the weights and the degen- 
eracy label. 

• class irrep: 

This class is the data structure of irreps that is needed to compute their 
tensor products. Its constructor takes as arguments the Lie algebra of type 
Liealg. liealgebra, the highest weight of type int list, and a third in- 
put parameter of type b * ((c * d) * e) where 
b = '_a list, 
c = (int * ket) list, 

d = ((int * int) * int Alg. lcrs_vector) list, 
e = ((int * int) * Alg. Icrs.f ield) list. 

In the case of non-degenerate or adjoint irreps consistent lowering normaliza- 



Note that in general degeneracies are not sufficient to single out adjoint represen- 
tations. 
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tions and scalar products can be computed as discussed in Section 2.3. The 
weight systems are obtained from the subprogram Liealg, cf. Section 4.2. 
In these cases the third input parameter can be chosen to be empty_input. 
On the other hand, more comphcated irreps first need to be constructed 
in a simpler tensor product. Subsequently, one reads off consistent lowering 
normalizations and scalar products which finally are to be used as an input 
parameter upon initializing an object of type irrep. In this case, the sec- 
ond argument specifying the highest weight is redundant and is ignored 
Here, type b becomes (int * lcrs_ketket) list, i.e. this part is the list 
of pairs of an integer enumerating all weights and the corresponding states 
in the tensor product representation of type lcrs_ketket (to be defined 
below). Type c is meant to model the list of pairs of the enumerating index 
and the ket (weights plus degeneracy index). Next, data of type d encodes 
the lowering normalizations: For every list entry, the second integer gives 
the simple root by which we lower the state specified by the first integer; 
the object of type int Alg. lcrs_vector gives the resulting linear combi- 
nations in terms of the unique state labels. Finally, type e encodes the scalar 
products: For every list entry, the two integers specify the states of which 
the scalar product is taken; the result is of type Alg. lcrs_f ield. 

Its methods are: 

get_la: -> Liealg. liealgebra: Returns the Lie algebra the irrep belongs 
to. 

get_hw: -> int list: Returns the highest weight of the irrep. 
get_dini: -> int: Returns the dimension of the irrep. 
is_adj oint : -> bool: Yields true if the irrep is adjoint and false otherwise. 
get_list: -> (int * ket) list: Returns the list that contains the unique 
labeling of the kets. 

print_list: -> (int * string) list: Like get_list but expresses kets 
as string. 

get_ket: int -> ket: Returns the ket specified by the integer. 
get_scalar_product : int -> int -> lcrs_field: 

Returns the scalar products of the states specified by the two integers. If the 
object of irrep has been initialized by means of a non-trivial third argu- 
ment resulting from another tensor product the scalar products are read off 
from this parameter. Otherwise, scalar products are computed as discussed 
in Section 2.3. 

lower: int -> int -> int lcrs_vector: Returns the linear combination 
that results from lowering the state specified by the second integer with the 
simple root denoted by the first integer. Again, lowering normalizations are 
read off from a non-trivial third input parameter of the constructor or fol- 
low from the rules for non-degenerate or adjoint representations discussed 
in Section 2.3 if the third input parameter is trivial. 



It suffices to put in an empty list [] . 
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class ['a] lcrs_2vectors: 

This class is an inheritor of the class [^a] lcrs_vector where 'a is taken 
to be ' a * 'a. Objects of this class can be viewed as elements of the tensor 
product of two ' a vector spaces. 

class lcrs_ketket: 

This class models states in the tensor product representation of two irreps. 
It is inherited from the class [int] lcrs_2vectors. Its constructor takes 
one argument of type (Alg. lcrs_f ield * (int * int)) list which is 
meant to be a linear combination over lcrs_f ield of basis states labeled by 
two integers which specify the respective states in the two irreps the tensor 
product is taken of. 

We have the following methods (in addition to the inherited ones) : 
add_lcrs_ketket : lcrs_ketket -> lcrs_ketket: 

weight : irrep -> irrep -> int list: Specifying the two irreps the ten- 
sor product is taken of, this method yields the weight of the state. Note that 
in order to save RAM an object of lcrs_ketket does not know which irrep 
it belongs to which is why such information must be given here (and in what 
follows) as an input parameter. 

print_kets: irrep -> irrep -> (string * (string * string)) list 

print_to_string: irrep -> irrep -> string: 

lower: b -> b -> int -> lcrs_ketket: where b abbreviates 

b = (int -> int -> int Alg. lcrs_vector) which is precisely the type 

of method lower of class irrep. In fact, the first two arguments should be 

supplied by the data from method lower of irrep for the two irreps of the 

tensor product. The third input parameter specifies the simple root that 

acts as lowering operator. 

scalar_product : lcrs_ketket -> b -> b -> Alg. lcrs_f ield: where b 
abbreviates b = (int -> int -> Alg. lcrs_f ield) which is the type of 
method get_scalar_product of class irrep. Also here, the second and the 
third argument should be supplied by the method get_scalar_product of 
the two irreps of the tensor product. The first argument is another object 
of lcrs_ketket with which the scalar product is to be computed. 

class irrep_in_product_rep : 

This class models irreps that appear in the tensor product. Its constructor 
takes only one argument of type lcrs_ketket which corresponds to the 
highest weight of the irrep. 

We have the following methods: 

get_hw: -> lcrs_ketket: Returns the highest weight. 

get_list : -> lcrs_ketket list list: Returns the product states in the 

irrep grouped according to their level. 
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get_dim: -> int: Returns the dimension of the irrep. 
print : irrep -> irrep -> 

(string * (string * string)) list list list: 
print_f lattened: irrep -> irrep -> 

(string * (string * string)) list list: 
dominant _we ight s : irrep -> irrep -> lcrs_ketket list: Returns the 
states whose weights are dominant. 

descend: irrep -> irrep -> irrep_in_product_rep: After instantiat- 
ing an object of class irrep_in_product_rep the full irrep is constructed by 
means of the method descend. Again the two irreps of the tensor product 
must be supplied as input parameters. 

prepare: irrep -> irrep -> b * ((c * d) * e) : Reads out consistent 
lowering normalizations and scalar products for an irrep in the tensor prod- 
uct and prepares this data as an input parameter for the constructor of class 
irrep. It is by this mechanism that one can in principle instantiate objects 
of class irrep that model arbitrary irreps. Here, type abbreviations read: 
b = (int * lcrs_ketket) list, 
c = (int * ket) list, 

d = ((int * int) * int Alg. lcrs_vector) list, 
e = ((int * int) * Alg. Icrs.f ield) list. 

In words, variables of these types contain the association of labels to states in 
the tensor product representatioi]0, the association of labels to the weights, 
consistent lowering normalizations, as well as the scalar products of the cho- 
sen basis states, respectively. 



• class clebsch_gordan_decomposition: 

This class is the data structure for the Clebsch-Gordan decomposition. It 
contains as its main part an implementation of the decomposition algorithm 
from Section 2.4. The constructor takes two arguments of type irrep that 
specify the two factors of the tensor product. Upon initialization only the 
highest weight representation in the product representation is computed. 

Its methods are: 

get_irrepa: -> irrep: Returns the first factor in the tensor product. 
get_irrepb: -> irrep: Returns the second factor in the tensor product. 
get_irreps_in_product : -> irrep_in_product_rep list: Lists all ir- 
reps that appear in the tensor product. 

basis_product : int list -> lcrs_ketket list: Computes a (non-ortho- 
gonal) basis for the tensor space for a given weight. 

decompose: -> unit: Completely decomposes the tensor product by con- 
structing states explicitly. 

check_dims : -> bool: Yields true if dimensions of the tensor product and 



' This bit of information is not needed in the constructor for class irrep which is 
why there this parameter is polymorphic. 
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its decomposition match, otherwise retm^ns false. 

result : -> string: Returns a string output of the decomposition. 

• class irrep_in_tp: 

This class has been designed in view of multiple tensor products: Its objects 
(corresponding to irrcps) can be tensorizcd to give a certain irrep in the 
tensor product which becomes again an object of irrep_in_tp. Hence, an 
iteration of this process is possible. We use a certain tree structure (int 
tree) the internal variable of this class in order to circumvent the type 
definiteness of classes in O'Caml. The constructor of this class takes an ar- 
gument of type irrep. 

Its methods are: 

get_irrep: -> irrep: Returns the irrep the object belongs to. 
get_res: -> (int * int tree Alg. lcrs_vector) list: 

Expresses states in this irrep in terms of the states of the multiple tensor 
product: It returns a list of pairs where the first entry labels the state in 
the irrep while the second entry in the pair is a linear combination of the 
tree that encodes the basis states of the multiple tensor product, 
otimes: irrep_in_tp -> int -> irrep_in_tp: Tensorizes the irrep with 
another irrep and returns another object of type irrep_in_tp that belongs 
to the irrep specified by the second argument according to the order in which 
irreps in the tensor product are constructed (1 gives the highest weight ir- 
rep, 2 the next one. . . ). 

expand: int -> int tree Alg. lcrs_vector: Returns the state specified 
by the first argument in terms of a linear combination of trees, 
untree: (int * string Alg. lcrs_vector) list: Lists all states in the 
irrep as a pair of the state label and a linear combination of bases in the 
tensor product which are represented as strings. 

filter: int -> int list -> irrep_in_tp: Filters out only terms that 
include product bases where numbers from int list appear in the place 
that corresponds to the tensor product factor specified by the first argu- 
ment. 

chbasis: int -> (int * int Alg. lcrs_vector) list -> 

irrep_in_tp: 

Performs a basis change in the tensor product factor specified by the first 
argument. The linear transformation is given by the second argument by 
means of associating to every state a linear combination of states. 
is_sym: int -> int -> int: Tests the symmetry under exchange of fac- 
tors in the tensor product specified by the first two arguments. Yields 
1,-1,0 in the case of symmetry, antisymmetry, or indefinite symmetry, 
respectively. 

scale : Alg. lcrs_f ield -> irrep_in_tp: Scales all states by the first ar- 
gument. 

tensor: int -> int list -> Alg. lcrs_f ield: Gives for the state spec- 
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ificd by the first argument, the coefficient that belongs to a product basis 
state specified by the second argument. 

Furthermore, we included some functions that can be useful in discussing 
symmetry breaking in particle physics by means of a generalized Higgs mech- 
anism: Invariant higher- dimensional operators (i.e. singlets in the muhiple 
tensor products) may contain what we caU Higgs fields which means that they 
are supposed to acquire a vacuum expectation value (vev) in a certain direction 
in the weight spaces that correspond to the irreps the Higgs transform under. 
Upon vev insertion higher-dimensional operators yield an effective potential 
which one may want to compute. The stabihty groups of those vevs determine 
the subsymmetry to which the fuU symmetry is said to be spontaneously bro- 
ken. It then makes sense to express the effective potential in terms of irreps of 
the subsymmetry. Here, we do not aim at giving a full set of routines that are 
capable of doing such computations in general. We rather give some general 
functions that we found useful and present one additional function that was 
well suited to our problem though its extension to similar problems should be 
possible. 

• e_lower: irrep -> int -> int Alg. lcrs_vector 

-> int Alg. lcrs_vector: This function models the lowering operator that 

is associated to the simple root specified by the second argument of type 
int. Given an irrep and a simple root it can be viewed as an operator that 
acts on an object of type int Alg. lcrs_vector and returns another object 
of this kind. 

• comm: (int Alg. lcrs_vector -> int Alg. lcrs_vector) 
-> (int Alg. lcrs_vector -> int Alg. lcrs_vector) 

-> int Alg. lcrs_vector -> int Alg. lcrs_vector: 
This function gives the commutator of two lowering operators. Such a func- 
tion is useful because lowering operators that correspond to arbitrary roots 
are in general given by such multiple commutators. For instance, one may 
want to break to a maximal subalgebra that is built from the zeroth root. A 
vev state should then be invariant under the multiple commutator of low- 
ering operators that corresponds to the zeroth root. 

• scp: irrep -> int Alg. lcrs_vector -> int Alg. lcrs_vector 

-> Alg. lcrs_f ield: Computes the scalar product of two states that be- 
long to a given irrep. 

• scalar_products : irrep -> int -> int 

-> int Alg. lcrs_f ield array array: Prints scalar products of weight 

states in irrep. The first integer argument specifies weight in irrep while the 
second integer specifies how many further weights are to be selected. Use 
Aux. print. array to nicely display result. 
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• gram: irrep -> int Alg. lcrs_vector list 

-> int Alg. lcrs_vector list -> int Alg. lcrs_vector list: 

Performs a partial Gram-Schmidt procedure. Given a subset of orthogonal 
states in a basis (second argument), it subtracts for the remaining states in 
the basis (third argument) their projections on the orthogonal states. Then, 
the original set of orthogonal states and the complement are orthogonal. 

• chbasis_list : int Alg. lcrs_vector list -> int 

-> (int * int Alg. lcrs_vector) list: The function chbasis_list re- 
turns the basis transformation rule in the format that is needed in the sec- 
ond argument of class irrep_in_tp, method chbasis. It computes how 
the generic weight states in the irrep, specified by an offset parameter of 
type int (second argument) and the length of the basis (first argument), 
are expressed in terms of the new basis vectors. 

• eff .couplings : irrep_in_tp -> irrep_in_tp 
-> (int list * string) array -> unit: 

This function is specialized to the case where one wants to compute the 
effective potential that results from an operator / multiple tensor product 
that consists of three identical irreps and another irrep that is said to ac- 
quire a vev which is left invariant by some subsymmetryf^ The output is 
the effective superpotential in terms of irreps of the subalgebra. The first 
argument is a singlet in the fourfold tensor product while the second one 
takes the singlet of the threefold tensor product (the fourfold product with- 
out the Higgs irrep). The latter is used as a means to recognize the terms 
that are invariant under the subsymmetry. The third argument is an array 
that contains the branching of the irrep with respect to the subsymmetry 
in that each array entry associates a label of type string to a subset spec- 
ified by the integer labels of weight states. The output consists of a linear 
combination of terms that are separately invariant under the subsymmetry. 
In this expression, every term is an abbreviation of the corresponding terms 
in the threefold tensor product (second argument). 



4.5 CGC^nh 



This part should be understood as a notebook file that is supplied with the 
computations one is interested in. Here, we include the three examples for 
testing purpose presented in Section 6. 



We were interested in the case of the tensor product 27 (g) 27 (8) 27 650 of 
that upon vev insertion yields an effective renormalizable superpotential. 
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5 Installation instructions 

5. 1 Toplevel mode 
Start O'Caml toplevel via: 

ocaml nums.cma Aux.cmo Liealg.cmo Alg.cmo CGC.cmo 

It is convenient to write operations in a notebook file "CGC_nb.ml" which can 
also be loaded into the topleveQ 

#use "CGC.nb.ml"; ; 

Functions / objects in modules "Aux", "Liealg", "Alg", "CGC", and "CGC.nb" 
can be accessed as follows: 

open Aux; ; 



5.2 Compilation mode 

Specify notebook file (e.g. CGC_nb.ml) and output executable in Makefile 
and run make. 

make 

Then, run executable. 

Note that the generation of weight systems (Liealg. ml) can be regarded as 
a subprogram of its own. As here input parameters are limited to the speci- 
fication of the Lie algebra and the heighest weight of the irrep, a standalone 
executable can be useful. We therefore included an option in the Makefile 
that compiles Liealg. ml and Liealg_ex.ml the latter of which specifies the 
input interface. The make command reads 

make lie . {opt I bin} 

depending on whether or not native- compiling is used. As a default, lie. opt 
is built; otherwise the user should edit the Makefile correspondingly. 



^ Of course, this can also be done in the step before. 
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6 Test run description 



6. 1 Derivation of explicit lowering normalizations for arbitrary irreps 



Bearing in mind that we only know consistent lowering normalizations and 
scalar products of basis states for non-degenerate or adjoint representations, 
see Section 2.3, we developed means to derive also those for more complex ir- 
reps: As, obviously, any irrep can be constructed in the (multiple) tensor prod- 
uct of non-degenerate or adjoint representations, basis choices, scalar products 
and their lowering normalizations can be analyzed therein. As an example, we 
give the O'Caml toplevel code needed to create an object associated to the 
650 of Eq which we know has many degeneracies. 



(* Create objects for the (anti-) fundamental E6 irreps, the 27 and 27bar *) 
let e6_27 = new irrep E6 [1;0;0;0;0;0] empty_input ; ; 
let e6_27bar = new irrep E6 [0;0;0;0; 1;0] empty_input ; ; 

(* Create object of the tensor product decomposition, 27x27bar *) 
let eg = new clebsch_gordan_decomposition e6_27 e6_27bar; ; 
(* ... computes the first irrep in 27x27bar, the 650 *) 
cg#decompose; ; 

(* ... computes all remaining irreps, yields: 
Dimensions match. 

Clebsch-Gordan decomposition successfully done! 

E6: (1,0,0,0,0,0, )27 x (0,0,0,0, 1 ,0, )27 = 

(1,0,0,0,1,0,)650 

(0,0,0,0,0,1,)78 

(0,0,0,0,0,0,)1 

*) 

(* Read out the 650 *) 

let tmp = nth cg#get_irreps_in_product 0;; 
(* Analyze lower ings and scalar products *) 
let low_scp = tmp#prepare e6_27 e6_27bar; ; 

(* Create object associated to the 650 irrep using low_scp as input *) 
let e6_650 = new irrep E6 [] low_scp; ; 

(* This object can now be used in further (multiple) tensor products *) 



34 



6.2 248^248 of Ei 



By the very nature of the exphcit decomposition algorithm, dimension of ten- 
sor products are hmited by memory and CPU time demands. While the 
650 §>> 650 of Eq exceeds 4 GB RAM, the software successfully produces 
Clebsch-Gordan coefficients for 248 248 of Es- Input code in an O'Caml 
toplevel is as follows: 

(* Create object for the fundamental and adjoint 248 irrep *) 
let e8_248 = new irrep E8 [0;0;0;0;0;0; 1;0] empty. input ; ; 

(* Create object of the tensor product decomposition *) 
let eg = new clebsch_gordan_decomposition e8_248 e8_248;; 
(* . . . computes the first irrep in 248x248, the 27000... *) 

cg#decompose ; ; 

(* ... computes all remaining irreps *) 
(* yields 
Dimensions match. 

Clebsch-Gordan decomposition successfully done! 

E8: (0,0,0,0,0,0,1,0,)248 x (0,0,0,0,0,0,1,0,)248 = 

(0,0,0,0,0,0,2,0,)27000 

(0,0,0,0,0,1,0 , 30380 

(1,0,0,0,0,0,0,0, )3875 

(0,0,0,0,0,0, 1,0, )248 

(0,0,0,0,0,0,0,0,)! 

*) 

(* Get singlet contraction in tensor product *) 
let singlet = nth cg#get_irreps_in_product 4; ; 
singlet#print e8_248 e8_248; ; 

(* 

[[[("-1", ("(0,0,0,0,0,0,1,0,)!", "(0,0,0,0,0,0,-1,0,)!")); 

("!", ("(0,0,0,0,0,1,-1,0,)!", "(0,0,0,0,0,-1,1,0,)!")); 

("-!", ("(0,0,0,0,1,-1,0,0,)!", "(0,0, 0,0, -1,1,0,0, )!")); 
...]]] 

*) 

6.3 Multiple tensor product: 4 (g) 4 (g) 6 (g) 15 of S'[/(4) 

In what follows we present the ' Caml toplevel code needed to compute the 
4-fold tensor products of SU (4) irreps, 

4(8)4(8)6(8) 15 = (35) 
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where . . . denote further non-singlet irreps. Subsequently, we let the 15 acqiiire 
an SU{3) x t/(l)-invariant vev and compute the resulting effective potential. 

(* Creating objects for irreps the tensor product is taken of *) 
(* As all irreps are either non-degenerate or adjoint, consistent 
lowering normalizations and scalar products can be generically computed *) 
let su4_4 = new irrep (A 3) [1;0;0] empty_input ; ; 
let su4_6 = new irrep (A 3) [0;1;0] empty_input ; ; 
let su4_15 = new irrep (A 3) [1;0;1] empty_input ; ; 

(* Creating objects for irreps that can be tensorized *) 
let t4 = new irrep_in_tp su4_4; ; 
let t6 = new irrep_in_tp su4_6;; 
let tl5 = new irrep_in_tp su4_15;; 

(* Define SU(3)xU(l)-singlet in 15 of SU(4) *) 

let sing_su4_15 = ((new lcrs_vector [(unity_f, 7); (number (-2) 1 1,8); 

(number 3 11, 9)])#scale (number 1 6 6) )#simplif y ; ; 
(* Check singlet, by descending with simple roots 1 and 2 *) 
(e_lower su4_15 1 sing_su4_15)#print ; ; 
(e_lower su4_15 2 sing_su4_15)#print ; ; 
(* In both cases, yields [] , i.e. *) 

(* Further states of weight zero *) 

let bsl5 = [new lcrs_vector [(unity_f, 7+1)]; 

(new lcrs_vector [(unity_f ,7+2)] )#scale (number 1 1 3)];; 
(* Check linear independence of states in the zero weight space *) 
linearly_dependent (sing_su4_15 : : bsl5) ; ; 
(* Yields false, i.e. vectors are linearly independent *) 
(* Using Gram-Schmidt procedure, choose basis in zero weight space 
consisting of the SU(3)xU(l) singlet and two states orthogonal to this 
singlet *) 

let zero_ws_su4_15 = gram su4_15 [sing_su4_15] bsl5;; 
let basis_su4_15 = [sing_su4_15] @ zero_ws_su4_15; ; 

(* Compute transformation of the generic basis to the one in basis_su4_15 *) 
let trafo = chbasis_list basis_su4_15 7; ; 

(* Compute all singlets in 4x4x6x15 *) 

(* t4#otimes t4 1 yields the first irrep in the tensor product 4x4 *) 
let ttl = ((t4#otimes t4 l)#otimes t6 2)#otimes tl5 7;; 

(* ... evaluation of the tensor products *) 

let tt2 = ((t4#otimes t4 2)#otimes t6 2)#otimes tl5 7;; 

(* ... evaluation of the tensor products *) 
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(* Check that tensor product yield singlets *) 
ttl#get_res ; ; 
tt2#get_res ; ; 

(* Yields [(1, <obj>)] => OK *) 

(* Check (anti) symmetry under exchange of first two irreps *) 

ttl#is_sym 12;; 

(* Yields 1, i.e. symmetric *) 

tt2#is_sym 12;; 

(* Yields -1, i.e. antisymmetric *) 



(* Filter out zero weight components 
4-fold tensor product *) 
(* Change zero-weight basis in 15 *) 
let ttl' = ttl#filter 4 [7;8;9];; 
let ttl" = ttl'#chbasis 4 trafo;; 
let tt2' = tt2#filter 4 [7;8;9];; 
let tt2'' = tt2'#chbasis 4 trafo;; 
(* Filter out vev component denoted 
let ttl'" = ttl"#filter 4 [-1];; 
let tt2'" = tt2"#filter 4 [-1];; 



(labeled 7,8,9) in the 15 in the 



1 *) 



(* Scale to nice overall normalization *) 
let ttl"" = ttl"'#scale (number 1 2 3);; 
let ttl"" = ttl"'#scale (number 1 1 3);; 



(* Print result *) 

(snd (hd ttl" "#untree))#print; ; 

(* Result: 

[("-1", "(((4,3),1),-1)"); ("-1", "(((3,4),1),-1)"); ("1", " (((4,2) ,2) ,-1)") ; 

("1", "(((2,4),2),-l)"); ("-1", "(((4,1),4),-1)"); ("-1", " (((1,4) ,4) ,-1) ")] *) 
(snd (hd tt2" "#untree))#print; ; 

(* Result: 

[("1", "(((1,3),5),-1)"); ("-1", "(((3,1),5),-1)"); ("-1", "(((1,2),6),-1)"); 
("1", "(((2,1),6),-1)"); ("1", "(((3,4),1),-1)"); ("-1", "(((4,3),1),-1)"); 
("-1", "(((2,4),2),-l)"); ("1", "(((4,2),2),-l)"); ("1", ••(((1,4),4),-1)"); 
("-1", "(((4,1),4),-1)"); ("-1", "(((2,3),3),-l)"); ("1", " (((3,2) ,3) ,-1) ")] *) 



(* Weight labels *) 
su4_4#print_list ; ; 
(* yields: 

[(1, "(l,0,0,)l"); (2, "(-l,l,0,)l"); (3, "(0,-1, 1,)1"); (4, "(0,0,-l,)l")] *) 



su4_6#print_list ; ; 
(* yields: 

[(1, "(0,1,0,)1"); (2, "(1,-1, 1,)1"); (3, " (1 ,0,-1 , ) 1") ; 
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(4, "(-1,0,1,)!"); (5, "(-1,1,-1,)1"); (6, "(0,-1,0,)!")] *) 
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A Appendices 



A.l Dynkin Diagrams 



Our conventions for Dynkin diagrams are shown in Fig. A.l. 



A. 2 Cartan Matrices 



We follow the conventions by Georgi [2] and Slansky [3] for the Cartan matrix, 
namely Aji = 2a^ ■ a^/la^Y which is the transpose of Cornwell [4]. For the 
Dynkin diagrams we follow the conventions of [2] concerning the normalization 
and numbering of the weights, which has n — 1 longer simple roots for S0{2n + 
1) and n—1 shorter roots for Sp{2n). For F4 the first two roots on the left 
of the Dynkin diagrams are the shorter ones, for G2 the first root on the left 
is the shorter, too. This is contrary to [4,3], where the shorter roots are on 
the right hand side. Therefore, the Cartan matrices for the infinite series and 
Eqj^s given here agree with [3] and are the transpose for F4, G2, while the 
opposite is true for [4]. Up to the definition of the Cartan matrix (which is 
important for the descending from the highest weight) and the refiection of 
the Dynkin diagrams for F4 and G2, we follow the conventions of Cornwell [4]. 
Our Cartan matrices are given by: 
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A. 3 Positive roots 



The positive roots are needed inside the Freudenthal formula for the calcu- 
lation of the degeneracy of the states in a representation. In this section we 
denote as usual the simple roots by a\ i — 1, . . . , rank(0). 



SU {n + 1) has the + 1) positive roots: 

k 

a" with j,k — 1,2, ... ,n; j < k 

u=j 

S0{2n + 1) has the positive roots: 

, j = 1,2, 
, j,k — 1,2, ... , 
, j,k ^1,2,... ,n j <k 



■IU=J 

-1 

=j 

-1 

■'U=j 



ES«" + 2E::=ik«« , j,k^l,2,...,n j <k 



Sp{2n) has the positive roots: 



j,k = l,2,...,n; j <k 



Ea«" + 2E;:=fea" + a^,j,k = l,2,...,n-l;j<k 



2e::4«"+«" 



, j = l,2,...,n- 1 
, j = l,2,...,n- 1 



• SO{2n) has the n(n — 1) positive roots: 
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n-2; j <k; 


(n > 3) 






+ «" 




1,2,... 


, n 


- 2 ; (n > 3) 












1,2,... 


, n 


- 2 ; (n > 3) 




E"~^ a" 






J = 


1,2,... 


, n 


- 2 ; (n > 3) 










j = 


1,2,... 


, n 


- 2 ; (n > 3) 





(A.9) 
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Eg has 36 positive roots (coefficients of tfie simple roots) : 



(100000), (011100), (mill) (010000), (ooiiio), (012101) 

(001000) , (001101), (112101) (000100), (111100), (012111) 
(000010), (011110), (112111) (000001), (001111), (012211) 
(110000), (111110), (112211) (011000), (011001), (122101) 

(001100), (111001), (122111) (000110), (011101), (122211) 

(001001) , (111101), (123211) (111000), (011111), (123212) 



E7 has the 63 positive roots (coefficients of tfie simple roots) : 



(1000000), 


(0001110), 


(0010000), 


(0111110), 


(0000100), 


(1111110), 


(0000001), 


(1110001), 


(0110000), 


(1111001), 


(0010001), 


(1111101), 


(0001100), 


(1111111), 


(0011001), 


(1121001), 


(0000110), 


(1121101), 


(0111100), 


(1121111), 


(1111100), 


(1122101), 



(0122111) 
(0122211) 
(1221001) 
(1221111) 
(1222111) 
(1232101) 
(1232211) 
(1232102) 
(1232212) 
(1243212) 
(2343212) 



(0100000), 
(0001000), 
(0000010), 
(1100000), 
(0011000), 
(1110000), 
(0111000), 
(1111000), 

(0011100) , 

(0011101) , 



(0011110) 
(0011111) 
(0110001) 
(0111001) 
(0111101) 
(0111111) 
(0121001) 
(0121101) 
(0121111) 
(0122101) 



(1122111) 
(1122211) 
(1221101) 
(1222101) 
(1222211) 
(1232111) 
(1233211) 
(1232112) 
(1233212) 
(1343212) 



£^8 has the 120 positive roots (coefficients of the simple roots) 



(10000000), 


(01111101), 


(12332101) 


(01000000), 


(11111101), 


(12332111) 


(00100000), 


(01111111), 


(12332211) 


(00010000), 


(11111111), 


(12333211) 


(00001000), 


(01210001), 


(12321002) 


(00000100), 


(11210001), 


(12321102) 


(00000010), 


(01211001), 


(12321112) 


(00000001), 


(11211001), 


(12322102) 


(11000000), 


(01211101), 


(12322112) 


(01100000), 


(11211101), 


(12322212) 


(00110000), 


(01211111), 


(12332102) 


(00100001), 


(11211111), 


(12332112) 


(11100000), 


(01221001), 


(12332212) 


(00011000), 


(11221001), 


(12333212) 


(01110000), 


(01221101), 


(12432102) 


(00110001), 


(11221101), 


(12432112) 


(11110000), 


(01221111), 


(12432212) 


(00001100), 


(11221111), 


(12433212) 


(00111000), 


(01222101), 


(12443212) 


(01111000), 


(11222101), 


(13432102) 


(00111001), 


(01222111), 


(13432112) 


(11111000), 


(11222111), 


(13432212) 


(00000110), 


(01222211), 


(13433212) 


(00011100), 


(11222211), 


(13443212) 


(00111100), 


(12210001), 


(13543212) 


(01111100), 


(12211001), 


(13543213) 


(00111101), 


(12211101), 


(23432102) 


(11111100), 


(12211111), 


(23432112) 


(00001110), 


(12221001), 


(23432212) 


(00011110), 


(12221101), 


(23433212) 


(00111110), 


(12221111), 


(23443212) 


(01111110), 


(12222101), 


(23543212) 


(00111111), 


(12222111), 


(23543213) 


(11111110), 


(12222211), 


(24543212) 


(01100001), 


(12321001), 


(24543213) 


(11100001), 


(12321101), 


(24643213) 


(01110001), 


(12321111), 


(24653213) 


(11110001), 


(12322101), 


(24654213) 


(01111001), 


(12322111), 


(24654313) 


(11111001), 


(12322211), 


(24654323) 
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Fi has the 24 positive roots (coefficients of the simple roots) 



(1000), (0011), (0121), (1222) 

(0100), (1110), (1121), (1231) 

(0010), (0111), (1220), (1232) 

(0001), (0120), (0122), (1242) 

(1100), (1111), (1122), (1342) 

(0110), (1120), (1221), (2342) 

G2 has the 6 positive roots (coefficients of the simple roots) : 

(10), (01), (11), (12), (13), (23) 
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